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ABSTRACT 




This report is Part H in a series of four reports on the statistical 
analysis of the Hawkins and Southworth random sample of 285 sporadic photo- 
graphic meteors. The four parts of the analysis comprise the combinations 
from two alternative formulations for meteoroid mass and two alternative for- 
mulations for data weighting. Parts I and II have the same weighting as a function 
of air-entry velocity and other parameters. But the surprisingly large disparity 
between the results with Opik's luminous efficiency suggests the possibility of 
considerable bias from weighting inversely with the 1. 5-power of velocity. 
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DEFINITION OF SYMBOLS 


Definition 

Meteoroid mass in grams in space 

Meteoroid velocity before deceleration in the atmosphere 

Celestial latitude of corrected radiant in degrees 

Meteoroid heliocentric orbital eccentricity 

Elongation of the true radiant from the apex of the earth's 
way in degrees 

Terrestrial weighting function in consideration of the population 
of meteoroids that actually encounter the earth 

Spatial weighting function in consideration of the population of 
meteoroids in space that may or may not encounter the earth 

Logarithm to base ten 

Representation of log m , log v, |0 e |> e > and X, respectively. 

Sample weighted mean, and weighted standard deviation for 
the variable X. where i = 1 5 

Sample weighted correlation coefficient between X. and X. 
for i, j = 1, .... 5 1 * 

The determinant of the sample weighted correlation coefficients 
r for i, j = 1, .... 5 

Cofactor of the element r„ in the determinant R 

ij 

Standard error of the estimate of X 4 from the equation of the 
regression plane 



DEFINITION OF SYMBOLS (Concluded) 


Symbol 


Definition 


r 1*2345 Multiple correlation coefficient for in relation to X 2 , . . . , X 5 


r 


ij*klm 


x 

□ 


M 

P 



F> 


F 

mv 


F 2 
mv 


F 

mv 


3/2 


Z 


Partial correlation coefficient between X. and X. with X, , X., 

1 J k 1 

and X held fixed, where i, j, k, 1 , m = 1 , . . . , 5 
m 

Data point for log m not less than the weighted median log m 
Data point for log m less than the weighted median log m 
Meteor absolute photographic magnitude 


Mean number of sporadic and stream meteors per second per 
square meter of level surface with absolute photographic 
magnitude equal to or less than M^ 

Mean number of sporadic and stream meteoroids per second 
per square meter of level surface with mass equal to or 
greater than m grams 

Mean number of sporadic and stream meteoroids per second 
per square meter of level surface with momentum equal to or 
greater than mv gram kilometers per second 

Mean number of sporadic and stream meteoroids per second 
per square meter of level surface with kinetic energy equal to 
or greater than mv 2 gram kilometers 2 second -2 

Mean number of sporadic and stream meteoroids per second 
per square meter of level surface with the geometric mean of 
momentum and kinetic energy equal to or greater than mv 3 / 2 
gram kilometers 3 / 2 second -3 / 2 

Zenith angle to meteor radiant in radians 
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STATISTICAL ANALYSIS OF PHOTOGRAPHIC METEOR DATA, PART H: 
VERNIANI'S LUMINOUS EFFICIENCY AND 
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SUMMARY 


This report is Part H in a series of four reports on the statistical 
analysis of the Hawkins and Southworth random sample of 285 sporadic photo- 
graphic meteors. The four parts of the analysis comprise the combinations 
from two alternative formulations for meteoroid mass and two alternative for- 
mulations for data weighting. Parts I and H have the same weighting as a function 
of air-entry velocity and other parameters. But the surprisingly large disparity 
between the results with Opik's luminous efficiency suggests the possibility of 
considerable bias from weighting inversely with the 1. 5-power of velocity. 


INTRODUCTION 


Justification and Purpose 

The technology of meteoroids and of their interaction with fields and with 
natural and artificial bodies in space is of considerable scientific and engineering 
interest. Each of the several sources of data continues to be of difficult inter- 
pretation because of indirection, extrapolation, and bias in random samples 
caused by physical selectivity. The difficulty of using a sample of photographic 
meteor data to infer the flux distributions of meteoroids as functions of dynamic 
parameters is that the results are sensitive both to the luminous efficiency 
formulation as a function of velocity and to the manner of weighting of the data. 

A study of the statistical consequences of such different alternatives should 
alleviate the decision problem. This report, with Part I [ i ] of this series , will 
present in more detail some results that were mentioned in two recent papers 
[2, 3]. 


Method and Notation 


A multivariable statistical analysis is made with weighting functions of 
meteor height, velocity, celestial latitude, zenith angle, and earth-encounter 
probability. The sample is then equally divided with respect to an intermediate 
mass value, and weighted cumulative distributions for several parameters are 
plotted for the two gradations with respect to mass, jointly or separately. The 
analysis is repeated without the weighting with respect to the earth-encounter 
probability. All weighting factors are adjusted so that the sum for the sample 
is equal to the sample size. All flux values are cumulative with respect to the 
indicated parameter (e. g. , mass, momentum, etc. ) and are in numbers per 
second per square meter of level surface. All logarithms are for base ten. 
Statistical notation is according to Hoel [4], 


Scope 

The same data sample and the same weighting as in Part I [1] of this 
series are used in the present analysis; i. e. , Hawkins and Southworth's [5, 6, ] 
random sample of 285 sporadic photographic meteors described in Part I [1], 

But instead of calculating the meteor mass values as in Part I [1] by using 
Opik's [7] physical theory of meteors, Hawkins and Southworth's 1 6] tabulated 
values are used. These mass values were computed by Hawkins and Southworth 
[6] using Hawkins [8] "short trail method," which is said to presuppose 
Verniani's [9] meteor luminous efficiency directly proportional to velocity. 

The weighting function for the present analysis was developed in Part I 
[ 1 ] . After weighting inversely with the square of meteor height at maximum 
brilliance, inversely with the 3/2 power of the air-entry velocity, inversely with 
the apparent fraction of the circle of celestial latitude through the meteor radiant, 
and inversely with Opik's [10] earth-encounter probability, it was found that a 
weighting also inversely with exp (0. 18Z) maximized symmetry with respect to 
the ecliptic plane (where Z is the zenith of the meteor radiant in radians). Be- 
cause the meteors with radiants more than 42° below the ecliptic were obscured 
by the horizon, the meteors with radiants more than 42° above the ecliptic were 
given double weight for arithmetic considerations of celestial latitude. This 
"spatial" weighting factor f g is used in consideration of a population of meteoroids 
in space regardless of whether or not they may encounter the earth. The "ter- 
restrial" weighting factor f t is similar to f g except that f t does not involve Opik's 
[10] probability that a meteoroid with given orbital parameters will encounter 
the earth during one revolution of the particle. The meteor data and values of 
the weighting functions were tabulated in Part I [ 1 ] . 
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DISCUSSION OF RESULTS 


The results with the set of mass values used in the present analysis are 
quite different from those reported in Part I [1] with the mass values calculated 
from Opik's [7] physical theory. The results of the multivariable statistical 
analysis are tabulated in Appendix I. For example, with "terrestrial" weighting, 
r^ (direct correlation between the logarithms of mass and velocity) changed 
from 0. 01 to -0. 69, while the corresponding partial correlation r^. 345 changed 
from 0. 10 to 0. 41; r M (direct correlation between log mass and eccentricity) 
changed from 0. 11 to -0. 45, while the corresponding partial correlation r 14 . 2 35 
changed only from -0. 14 to -0. 23; and r 15 (direct correlation between log mass 
and elongation from the apex of the earth’s way) changed from 0. 11 to 0. 54, 
while the corresponding partial correlation r 15 . 234 changed only from -0. 02 to 
-0. 07. The indicated changes in r 14 and r 15 appear to reflect the changes in r^ 
and r 12 . 345 through r 24 = 0. 81 and r 25 = - 0. 64. 


With "spatial" weighting, the median value log m = -1. 88 divides the 
sample into two parts with a lower weighted mean log m = -2. 35 and an upper 
weighted mean log m = -1. 38. With "terrestrial" weighting the corresponding 
median, lower mean, and upper mean for log m are -1. 46, -1. 98, and -1. 17, 
respectively. Figures 1 through 5 each show cumulative distributions of para- 
meters over the separate mass-subsets. 

Cumulative distributions of eccentricity are shown in Figures 1 and 2 
for "spatial” and "terrestrial" weighting, respectively. The separation of the 
"high" mass cumulative distribution from that for the meteors of "low" mass 
depends more on the ordinary correlation r 14 than on the corresponding partial 
correlation r I4 . 235; e. g. , Figures 1 and 2 show that meteors of low mass tend 
to have high eccentricity in agreement with the negative correlations r 14 = -0. 47 
and -0. 45, respectively, from Appendix I. 

Figures 3 and 4 show for the cumulative distribution of arithmetic 
celestial latitude of meteor radiant the same as Figures 1 and 2, respectively, 
showed for the cumulative distribution of eccentricity. The barely separable 
plots in Figures 3 and 4 reflect the numerically low correlations r 13 = 0. 10 and 
0. 04, respectively, from Appendix L But the widely separated "terrestrially" 
weighted cumulative distributions of velocity in Figure 5 reflect the strongly 
negative correlation r 12 = -0. 69 from Appendix I. 
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► 


Figures 6 through 10 show weighted whole -sample log cumulative distri- 
butions with respect to the logarithms of dynamic parameters such as mass, 
momentum, etc. In each case the plots seem to be approximately linear over 
the large-mass half of the sample weight. Presumably, the smaller masses are 
not adequately represented, and therefore should be ignored in Figures 6 through 
10 . 


Figure 6 shows that the logarithm of the "spatially" weighted cumulative 
distribution of log mass has a unit negative slope with respect to log mass for 
the present set of mass values, just as was found in Part I [1] with the other 
set. Also the corresponding slope with "terrestrial" weighting in Figure 7 has 
the same value (-1. 34) as was found in Part I [1] , with the other set of mass 
values, and as was found previously by Hawkins and Upton [11] with a somewhat 
smaller sample and different weighting. 

Figures 8 through 10 show that the slope (-1. 34) in Figure 7 is invariant 
with respect to replacing log mass with the logarithms of momentum, kinetic 
energy, and the geometric mean of momentum and kinetic energy, respectively. 
The derivation in Appendix II indicates that this is the result that should be found 
if mass and velocity were statistically independent. But the correlation is not 
numerically small in the present case; and different slopes were found in Part I 
[1] with numerically smaller correlation. No explanation for this effect has 
been found. 

Any ordinate in Figures 7 through 10 is converted into log cumulative 
mean total flux by subtracting -15. 18 (the logarithm of the area-time exposure 
product in square meter seconds) and adding 0. 08 (the logarithm of the factor 
by which mean total flux exceeds mean sporadic flux) ; see Part I [ 1]. The 
results are shown in Figure 11. 


CONCLUSIONS 


The results of the present analysis differ more. "widely from those in 
Part I [ 1] with the same weighting than had been expected. This is caused by 
the interacting roles of mass and velocity, and possibly to a bias with respect 
to velocity that may not be appropriately reduced by the weighting function that 
is inversely proportional to the 3/2 power of velocity. The next run of this 
analysis, with the Hawkins and Upton's [11] weighting as a function of magnitude 
above the limit of the photographic plate , instead of weighting as a function of 
velocity, is now expected to give more convincing results. 
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FIGURE 1. SPATIALLY WEIGHTED DISTRIBUTIONS OF ECCENTRICITY 
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FIGURE 2. TERRESTRIALLY WEIGHTED DISTRIBUTIONS OF ECCENTRICITY 
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Sample Cumulative Distribution Separately for 'High* and 'Low 
Mass Regimes with Spatial Weighting f s 
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FIGURE 3. SPATIALLY WEIGHTED DISTRIBUTIONS OF CELESTIAL 
LATITUDE OF RADIANT 




Sample Cumulative Distribution Separately for 'High' and 'Low 
Mass Regimes with Terrestial Weighting ff 
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FIGURE 4. TERRESTRIALLY WEIGHTED DISTRIBUTIONS OF 
CELESTIAL LATITUDE OF RADIANT 



















Sample Cumulative Distribution Separately for 'High' and 'Low' 
Mass Regimes with Terrestrial Weighting f t 
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FIGURE 5. TERRESTRIALLY WEIGHTED DISTRIBUTION OF 
AIR-ENTRY VELOCITY 


















Logarithm of Sample Cumulative Distribution 
with Terrestrial Weighting 
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FIGURE 7. TERRESTRIALLY WEIGHTED DISTRIBUTION OF METEOROID MASS 
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FIGURE 8. TERRESTRIALLY WEIGHTED DISTRIBUTION OF 
METEOROID AIR-ENTRY MOMENTUM 
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FIGURE 10. TERRESTRIALLY WEIGHTED DISTRIBUTION OF THE GEOMETRIC 
MEAN OF METEOROID AIR-ENTRY MOMENTUM AND 
KINETIC ENERGY 
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tog F Mp * 0.537 Mp -13.81 
log F,. * -1.34 log m -14.94 
log F mv « -1.34 log (mv) -13.34 
log F mv z = -1.34 log (mv 2 ) -11.62 
log F mv s/2 * -1.34 (mv 3/z ) -12.38 
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PARTIAL CORRELATIONS 


Terrestrial Weighting f t 

0.405 for log m vs. log v 

0.195 for log m vs. |j3 e | 

-0.435 for log v vs. |£e| 

0.835 for log v vs. X 


Spatial Weighting f s 

0.135 for log m vs. |j9 e | 
-0.235 for log m vs. e 
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APPENDIX I. NUMERICAL RESULTS FROM THE MULTIVARIABLE 

STATISTICAL ANALYSIS. X 4 X 5 = log m, log v, 

1/3 I, e, AND X, RESPECTIVELY 


STATISTICAL 

PARAMETER 

WITH SPATIAL 
WEIGHTING f 

s 

WITH TERRESTRIAL 
WEIGHTING f 

Xi 

-1. 8652 

-1. 5691 

x 2 

1. 4240 

1. 2870 

X 3 

28. 9383 

25. 8878 

X 4 

0. 7040 

0. 5505 

X 5 

83. 1402 

97. 1172 

si 

0. 6285 

0. 5433 

s 2 

0. 1793 

0. 1684 

s 3 

18. 9195 

22. 2727 

s 4 

0. 2214 

0. 2418 

s 5 

21. 2353 

24. 4780 

r ll 

1. 0000 

1. 0000 

r 12 

-0. 7101 

-0. 6879 

r 13 

0. 1040 

0. 0380 

r 14 

-0. 4714 

-0. 4529 

r l5 

0. 6828 

0. 5429 

r 22 

1. 0000 

1. 0000 

r 23 

-0. 1242 

0. 0465 

r 24 

0. 8324 

0. 8108 

r 25 

-0. 8258 

-0. 6369 

r 33 

1. 0000 

1. 0000 







STATISTICAL 

WITH SPATIAL 

WITH TERRESTRIAL 

PARAMETER 

WEIGHTING f 

s 

WEIGHTING f 

L 

r 34 

-0. 2718 

-0. 1343 

*35 

0. 1196 

-0. 0606 

*44 

1. 0000 

1. 0000 

r 45 

-0. 4496 

-0. 1649 

*55 

1. 0000 

1. 0000 

R 

0. 0132 

0. 0313 

R 11 

0. 0300 

0. 0658 

Rl2 

-0. 0361 

-0. 0763 

R 13 

-0. 0032 

-0. 0100 

R 14 

0. 0155 

0. 0323 

R 15 

0. 0027 

0. 0070 

R22 

0. 3690 

0. 5399 

^23 

0. 0392 

0. 0637 

R24 

-0. 2166 

-0. 3726 

**25 

-0. 1872 

-0. 2371 

R 33 

0. 0185 

0. 0396 

R34 

-0. 0270 

-0. 0483 

R 35 

-0. 0202 

-0. 0248 

R 44 

0. 1459 

0. 2976 

R 45 

0. 1059 

0. 1678 
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STATISTICAL 

PARAMETER 

WITH SPATIAL 
WEIGHTING f 

s 

WITH TERRESTRIAL 
WEIGHTING f 

R 55 

0. 1210 

0. 1493 

s 

0. 4173 

0. 3748 

e 

r 1* 2345 

0. 7478 

0. 7240 

r 12* 345 

0. 3426 

0. 4051 

r 13* 245 

0. 1353 

0. 1951 

r 14* 235 

-0. 2346 

-0. 2308 

r 15* 234 

-0. 0445 

-0. 0704 

r 23* 145 

-0. 4748 

-0. 4352 

r 24* 135 

0. 9332 

0. 9296 

r 25* 134 

0. 8870 

0. 8350 

r 34* 125 

0. 5209 

0. 4451 

r 35" 124 

0. 4282 

0. 3219 

r 45* 123 

-0. 7966 

-0. 7958 








APPENDIX II. SLOPE OF LOO-CUMULATIVE-FLUX VS LOG-MASS, 
LOG-MOMENTUM, ETC. FOR METEOROIDS ASSUMING 
STATISTICAL INDEPENDENCE OF MASS AND VELOCITY 


Let F ># the meteoroid cumulative-flux with respect to mass m, be 
F > = llA i / 2 , ( 1) 


where /? 6 and /3 2 are constants. A meteoroid with mass m not less than some 
limiting value m^ but otherwise random, has the following probability density 
function for m: 


f(m) 



( 2 ) 


A meteoroid with m > m^ that has a particular velocity v, but that is otherwise 
random, can satisfy the requirement 

mv Q > M (3) 


with probability P , 
c 


/ -n \^ 2 

P = f°° f(m)dm = (^— ) , I 

C Mv _n \ m L / 

where n is a constant; i. e. , n = 1 for momentum, 2 for kinetic energy, etc. 
Then, if v is randomly distributed in the interval v^ < v < v u and has some 
continuous probability density function f(v) , then the probability P that a 
meteoroid with m > but otherwise random will satisfy mv > M is found 
from 
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where k is constant is determined by the further condition that P = 1 when 
„ n . 

M = m v ; 1 . e. , 

Li Li 




(7) 


Therefore, because M in Equation (7), and m in Equation (1) have the same 
exponent, the slope of log-cumulative-flux vs log-mass should be invariant with 
respect to the substitution of momentum or energy for mass when mass and 
velocity are statistically independent. 
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